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Abstract. We have successfully ported under GNU Linux ATLAS 9, the widely 
used stellar atmosphere modeling code, as well as both the SYNTHE suite of pro- 
grams, its "companion" for spectral synthesis, and WIDTH, used to derive chem- 
ical abundances from equivalent widths of spectral lines. The porting has been 
realized by using the Intel Fortran Compiler. Our aim was to port the codes with 
the minimum possible amount of modifications: full compatibility with the VMS 
version has been maintained, along with all the codes functionalities. Dramatic im- 
provement in calculation speed with respect to the VMS version has been achieved. 
The full suite of codes is intended to be freely available to anyone. 
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1. Introduction 



Version 9 of the ATLAS code ("Kurucz ' 
1X97(1 [12^ is widely used for the produc- 
tion of LTE one-dimensional atmosphere 
models. Its fields of application extend from 
the modeling of stellar atmosphere for the 
purpose of stellar abundance analysis, to 
the computations of extended grids of col- 
ors for the production of synthetic CMDs, 
to the computation of large grids of syn- 
thetic spectra for stellar population synthe- 
sis. ATLAS is complemented by a suite of 
codes for spectrum synthesis (to which we 
shall refer collectively as SYNTHE) and by 
WIDTH, a code to derive chemical abun- 
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dances from the equivalent widths (EW) of 
observed lines. Since now, the original ver- 
sions of the codes were developed only for 
VMS machines, although some unix port- 
ings exist. Most of these portings, never- 
theless, depart significantly from the origi- 
nal code, and have tipically been adjusted 
to the specifical needs of the people in- 
volved in the porting. We undertook the 
porting under GNU Linux of ATLAS 9, 
WIDTH and SYNTHE packages with two 
main goals in mind: 

— To make the codes available under 
Linux, an operating system nowadays 
far more widespread than VMS, and 
running on much less expensive ma- 
chines; 

— To port the code introducing the small- 
est possible alterations in it, leaving 
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ATLAS 


SYNTHE 


VMS 

Pentium 4 
Pentium M 


478 (3.54) 
177 (1.31) 
122 (0.9) 


69 
10 
10 



Table 1. Execution time comparison for 
ATLAS 9, calculated on 135 iterations of a 
72 layers model (time per iteration within 
parentheses), and for the calculation of 
a 5 nm synthetic spectrum at resolution 
600000 with SYNTHE. Systems are: VMS: 
AlphaServer 800 5-500; Pentium 4: 1.9 GHz 
768 Mb RAM kernel 2.4.18-3 (Red Hat 
7.3) IPC 7.0; Pentium M: 1.6 GHz, 512 
Mb RAM kernel 2.4.22-1.2115.nptl (Fedora 
Core 1) IPC 8.0. All times in seconds. 



unchanged input/output formats and 
functionalities of the programs. This 
would allow a smooth integration with 
the already existing line databases, 
opacities and model grids. Moreover, 
this would allow to move effortlessly to 
Linux any new modification in the VMS 
code. 



2. Porting with Intel Fortran 
Compiler 

ATLAS 9 and the other codes in the suite 
have been developed along some 40 years. 
They are mostly written in Portran IV. 
Thus, a very good set of backward compat- 
ibility options were needed in the compiler, 
in order to minimize the needed changes. 
We choose to employ the Intel Portran 
Compiler (IPC), a P90-95 (with proprietary 
extensions) compiler whose Linux version 
is licensed free of charge for research ap- 
plication. Versions 7.0 and 8.0 have been 
used. The changes to the codes were mini- 
mal. Among the most important: 

— Some FORMAT statement syntaxes have 
been updated; 



Fig. 1. Plot of T against log{Tjioss) along 
the atmosphere of a star of T^ff = 5000 K, 
log g = 2.5, [Fe/H]=-0.5, as produced by 
the Linux ported code (solid line) and the 
original VMS code (triangles). 



The syntax of the OPEN commands have 
been modified to fit the new UNIX ex- 
ecution scripts; 

Some variables have been redefined to 
keep them in double precision through- 
out the various routines; 
Some big data tables (e.g. partition 
functions, integration matrices) were 
originally written in form of huge DATA 
blocks. IPC appears to have a limit to 
the size of DATA blocks, so we needed 
either to move them to files or to split 
them into smaller DATA blocks to avoid 
compilation problems. This lead to the 
need for ATLAS to read a further in- 
put file (PFIRON.DAT) with respect to 
the VMS version; 

Molecular input data in SYNTHE have 
been passed to ASCII format (from 
the VMS binary of the original). The 
program reading these files has been 
changed accordingly. 
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Fig. 2. Synthetic spectrum for the above mentioned star around the Mg b triplet. Open 
squares is VMS original code, solid line Linux code. The resolution of the calculation is 
600000, then a gaussian broadening (7 km/s FWHM) has been added to simulate the 
output of a spectrometer with resolution of about 43000. An example of line label: 518360 
= 518.360 nm; 1200 = Mg (12 atomic number) neutral (00); 2 = 2% residual intensity 
at line center (unbroadened spectrum); AZ is a code for the line \og{gf) source. 



3. Scripts and input data 

Original VMS execution scripts have ob- 
viously been rewritten as C- shell scripts, 
but the program sequence and inputs re- 
main unchanged (except for the need to 
read PFIRON.DAT). AU the ASCII inputs 
(models, line lists, Rosseland opacities) re- 
maining the same, we needed to "trans- 



late" the Opacity Distribution Functions 
(ODF), since they were written in VMS bi- 
nary form, unreadable under Linux. They 
have been copied in ASCII format and then 
rewritten in Linux binary (ODFs are bulky 
and the binary format saves significant disk 
space). We have p orted pro g ram v ersions 
for both the "old" iKiirucz I iITqqI ODFs 
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Fig. 3. Same as in fig. Q] but now compar- 
ing electron number densities. 



and the 
ODFs. 



"new" ICastelli fc Kunicz I (|2nn,'^ 



4. Performance 

The rcsuhs obtained with the ported ver- 
sion of the codes are, as can be seen in fig. 1 
through 3, indistinguishable from the ones 
produced by the VMS version. 

The porting leads to a dramatic in- 
crease in performance with respect to the 
VMS version. Examples of ATLAS and 
SYNTHE execution times (WIDTH is a 
very fast program, with execution times of 
about 1-2 seconds) are summarized in ta- 
ble Q Listed execution times are indicative 
since they also depend on disk throughtput. 
As can be seen, both atmosphere model- 
ing and spectral synthesis can now be per- 
formed in very short times on a mainstream 
laptop computer (the system labeled as 
"Pentium M" in tabled is actually a lap- 
top). This would allow any user to calculate 
directly any atmosphere he needs, without 
having to rely on precalculated grids. 

Moreover, together with the use of 
Linux supercomputer architectures such 
as Beowulf, this allows to use this pro- 



grams in the framework of large automated 
spectral analysis codes (see for instance 
iBonifacio fc CafFau1l2003(l conceived to an- 
alyze semi-automatically the high amount 
of data produced by high multiplexity 
spectrographs such as VLT-FLAMES. The 
high speed at which both the atmosphere 
modeling and the spectral synthesis steps 
can be performed, allows to produce large 
amounts of Monte Carlo tests to assess the 
error budget of the derived abundances. 

5. Availability 

The suite, as the original Kurucz one, is 
intended to be freely available to every- 
one interested in it. The full disk space re- 
quired is about 700 Mb, for a typical set 
of ODFs of the BIG type (for model com- 
putation), line lists and molecular data. 
For the time being, no webpage has been 
set to download the codes from, so inter- 
ested people are encouraged to contact us 
to obtain the codes and any information re- 
lated to them (email IsbordonOeso . org or 
bonif aciSts . astro . it ). 

Anyone using the codes for scien- 
tific publications is only requested to cite 
this work along with Kurucz (1993). 
Nevertheless, we plan to publish it on the 
Web as soon as possible, along with the 
proper documentation. 
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